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ABSTRACT 

We calculate the expected Mid-IR molecular hydrogen line emission from the first 
objects in the universe. Due to their low masses, the stellar feedback from massive 
stars is able to blowaway their gas content and collect it into a cooling shell where 
H2 rapidly forms and IR roto- vibrational (as for example the restframe 2.12 /im) lines 
carry away a large fraction (up to 10%) of the explosion energy. The fluxes from these 
sources are in the range 10~ 21 — 1CP 17 erg s _1 cm~ 2 . The highest number counts 
are expected in the 20 /mi band, where about 10 5 sources/deg 2 are predicted at the 
limiting flux of 3 x 10~ 18 erg s _1 cm~ 2 . Among the planned observational facilities, 
we find that the best detection perspectives are offered by NGST, which should be 
able to reveal about 200 first objects in one hour observation time at its limiting flux 
in the above band. Therefore, Mid-IR instruments appear to represent perfect tools 
to trace star formation and stellar feedback in the high (z ^ 5) redshift universe. 
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1 INTRODUCTION 

Detecting the first luminous objects in the universe will be 
the primary goal of several future space- and ground-based 
telescopes. The importance of these experiments consists in 
the fact that they could be able to test current cosmologi- 
cal scenarios, study the properties of these (supposedly very 
small) galaxies, and their effects on the surrounding envi- 
ronment (as for example reionization, heating and metal 
enrichment of the IGM). Due to the predicted low lumi- 
nosity of such objects, this task will be at the capability 
edge of even the most advanced and powerful instruments. 
Apart from the indirect probes of their effects, as for ex- 
ample the secondary anisotropics in the CMB left by reion- 
ization (Knox, Scoccimarro & Dodelson 1998; Gruzinov & 
Hu 1998; Bruscoli et al. 2000; Benson et al. 2000), some 
search strategies have already been suggested in the litera- 
ture. Marri & Ferrara (1998) and Marri, Ferrara & Pozzetti 
(2000) have suggested that Type II supernovae occurring in 
the first objects might outshine their parent galaxy by more 
than a hundred times and become visible by instruments 
like the Next Generation Space Telescope (NGST). Schnei- 
der et al. (2000) have investigated the possibility of detecting 
gravitational wave emission from high-redshift very massive 
objects with LISA. Oh (1999) proposes direct imaging of the 
ionized halos around primordial objects either via their free- 
free emission (possibly detectable with the Square Kilometer 
Array) or Balmer line emission again with NGST. 



One point that is particularly important when dealing 
with these low mass systems is that they are strongly af- 
fected by feedback mechanisms both of radiative and stel- 
lar type; these have been extensively investigated by Ciardi 
et al. (2000). Molecular hydrogen, being the only available 
coolant in a plasma of nearly primordial composition, is a 
key species in the feedback network as it regulates the col- 
lapse and star formation in these objects. Ferrara (1998, 
F98) pointed out that H2 is efficiently formed in cooling 
gas behind shocks produced during the blowaway (i.e. the 
complete ejection of the galactic gas due to SN explosions) 
process thought to occur in the first objects, with typical H2 
fractions fu 2 ~ 6 x 10~ 3 . We shall see that the conditions in 
these cooling blastwaves are such that a noticeable amount 
of the explosion energy is carried away by infrared (red- 
shifted into the Mid-IR [MIR] spectral region) H2 molecular 
lines, which therefore might provide us with a superb tool 
to detect and trace these very distant primordial galactic 
blocks. 

The use of molecular lines as diagnostics of moder- 
ate redshift (z ^ 3) galaxies has already been proposed by 
some authors (Frayer & Brown 1997; Blain et al. 2000); at 
higher redshifts pioneering calculations were carried on by 
Shchekinov & EnteP (1985); more recently Silk & Spaans 
(1997) concentrated on the CO and dust emission from HII 
regions inside larger galaxies. All these studies have empha- 
sized the power of the molecular line emission as a probe of 
distant sources. Based on previous calculations (F98), in this 
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paper we calculate MIR fluxes and number counts for these 
sources in various bands and assess if forthcoming observing 
facilities will be able to eventually unveil the beginning of 
the cosmic star formation era and study the feedback pro- 
cesses in the young universe. 



2 EMISSION MODEL 

Due to its symmetry, the H2 molecule has no electric dipole 
moment in the ground state. Therefore, the first detectable 
H2 emission lines are produced by quadrupole radiation and 
they are purely rotational. The line properties of interest to 
the present work, both for rotational and roto-vibrational 
lines, are listed in Table 1. There we see that the excitation 
temperatures, T ex , of these lines fall in the range 500-7000 
K. The temperature range spanned during the cooling of 
the post-shock IGM gas produced during the blowaway of 
low-mass primordial galaxies is (F98): 



Table 1: Relevant H2 Transition Lines 



300 K < T < 2.3M 6 2/5 (1 + z 



\18/5 



K 



(1) 



where M — 10 6 M6 Mq is the total mass of the galaxy.* 
The upper limit is set by the condition that the cooling 
time of the post-shock gas is shorter than the Hubble time. 
Strictly speaking the above upper limit holds until the IGM 
cooling is dominated by the inverse Compton radiative losses 
(z £ 6), but we use it down to z = 4, where it still represents 
a very good approximation for our purposes. Given the range 
in eq. 1, it is then conceivable that the above molecular lines 
are excited during the process in which a cold shell forms, 
thus producing potentially detectable radiation. The flux 
observed on the ground in a given line is derived as follows. 

Let Vik be the restframe frequency of photons emitted 
by the molecule during the transition between the energy 
levels k and i; then the line emissivity (erg s _1 cm~ 3 sr _1 ) 
is (Spitzer 1978): 

3»ik = ^v lk n k A k i, (2) 

where njt is the number density of molecules in the k level 
and Aki is the Einstein coefficient for spontaneous emis- 
sion. The typical H2 densities found are much lower than 
the critical one (~ 10 4 cm -3 ), hence we neglect collisional 
de-excitations. 

As the conditions for thermodynamic equilibrium are 
not satisfied, the population of the various levels must be 
obtained by solving the detailed balance equations (Spitzer 
1978). Using this approach, the number density of molecules 
in the vibrational level v, nj, of even and odd rotational 
levels decouple, and they are obtained by iteration using 
the following formula: 

-AE/kT 



n J+ 2(v) = n.j(v) gj 



(3) 



1 + 7J 

where gj is the statistical weight of level J, 7j is the ratio 
between the collisional excitation rate and the Einstein A 
coefficient, AE is the energy difference between the two lev- 
els. Obviously, in the limit of large 7j eq. 3 approaches the 
Boltzmann distribution. 

* We adopt a ACDM (cluster normalized) cosmology with f2j^f = 
0.35, fl A = 0.65, fl b = 0.04, h = 0.65. 
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1015 
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0-0S(3) 
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6.9 
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5.88 xlO- 8 


1-0S(1) 


2.12 


6953 


3.47 xl0~ 7 



(^Emission wavelength 

( 2 ) Excitation temperature (Combes & Pfenninger 1998; Timmer- 
mann et al. 1996) 

( 3 ) De-excitation Einstcn coefficient (Turner, Kirby-Docken & 
Dalgarno 1977) 



Given the post-shock temperatures found (eq. 1), it is 
only necessary to consider the two vibrational levels v = 0, 1. 
The relative number of molecules in these two levels depends 
on the total molecular hydrogen density, uh 2 ■ For a given 
redshift, the value of uh 2 (F98) is: 



n,H 2 (z) = p }h 2 (M, 



■PfH 2 {M,z)6 n H (z), 



(4) 



where fu 2 is the molecular fraction in the shell, and uh(z) 
is the IGM hydrogen density. We allow for a density en- 
hancement in the shell with respect to nn equal to <5; this 
is produced both by the shock compression and the possible 
occurrence of the explosion inside an overdense region of the 
universe (e.g. a cosmological filament). Because of the first 
effect, and as the expanding blastwave becomes rapidly ra- 
diative, 8 is given by the square of the shock Mach number, 
M, with respect to the ambient IGM. Following F98, it is 
easy to show that the density enhancement is given by: 



^4.8M 6 2/5 (l + z) 13/5 . 



(5) 



For simplicity, we have neglected the additional density in- 
crease produced by an explosion in a overdense region; thus, 
the results presented above should be interpreted as a con- 
servative lower limit to the detectability of the sources under 
investigation. Finally, we assume a branching ratio p = 0.75 
(0.25) for ortho (para) transitions which is valid under LTE 
conditions; however, non-equilibrium conditions might lead 
to slightly lower values of the ortho-to-para ratio (Chrysos- 
tomou et al. 1993; Rodriguez- Fernandez et al. 2000). 

With the emissivity given by eq. 2, the observed flux 
(erg s~ 



1 cm 2 ) is: 



.(1+*) 



(6) 



where v = fik/{^- + z) is the observed frequency, L Vik the 
luminosity and d L is the cosmological luminosity distance. 
For a flat universe (Oq = Qm + = 1), dz, is given by 



d L = {l + z) ( (1 + z') 
Jo 



dt 



dz> 



dz , 



(7) 



H (l + z)y/{l + Q M z)(l + z) 2 - Q A z(2 + z), (8) 



where H = lOO/i km s 1 Mpc 1 is the current Hubble 
constant. The corresponding luminosity is: 



(9) 
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where V = (4/3)ivRgdR a is the physical volume occupied by 
the H2 forming shell. The values of R s and dR s are obtained 
from the formulae in F98; the only difference consists in the 
assumption of a single burst of star formation, as opposed 
to the quiescent star formation prescription of F98. By as- 
suming a star formation efficiency of 10%, and one SN every 
100 Mq of stars formed, it is easy to show that the total 
energy of the explosion is equal to E — 1.3 x 10 53 M6 erg. 
Then, from F98 it follows: 

R s ~ 0.88A/ 6 1/5 (1 + Z y 11/5 Mpc, (10) 
dR s = ^ ^ (11) 

where we have assumed that all the swept IGM mass is in 
the cool shell of total mass M a , together with the galaxy 
interstellar medium. This result is somewhat dependent on 
the assumption made to derive the energy of the explosion, 
while it is not strongly dependent on the cosmology. 
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3 NUMBER COUNTS 

The number of objects whose observed flux is larger than 
the threshold value Fmin i s: 



N(> Fmin) = 



dz-^- 

dz 



dL 



dz 



dV 

dz 



dM 



dn dM 
'dM~dL 

dn 



M min 



dM' 



(12) 



where dn/dM is the comoving number density of halos with 
masses between M and M + dM (Press & Schechter 1974) 
and dV/dz is the comoving volume element per unit redshift: 



dV_ 

dz 



Aivcdf 
(1+^) 



(13) 



where dh and dt/dz are given in eqs. 7 and 8, respectively. 
The limits z m in and z ma x depend on the observational wave- 
length band; L ma x is the luminosity corresponding to the 
maximum mass value, M max , experiencing blowaway at red- 
shift z (Ciardi et al. 2000); M min = mm[M H , M(L min )}. 
Here Mh is the minimum halo mass in which stars can form 
at a given redshift, which is related to the minimum virial 
temperature of w 10 4 K below which atomic hydrogen cool- 
ing of the gas is suppressed and fragmentation into stars is 
inhibited (Haiman, Rees & Loeb 1997; Ciardi et al. 2000); 
M(L m i n ) is the mass of the halo producing the luminosity 
L m in corresponding to the observed flux F m in- Finally, the 
luminosity, L, is taken as the maximum one reached by the 
object (usually corresponding to temperatures equal to the 
excitation ones) that experiences a blowaway at redshift z, 
providing that it is emitted at a redshift > z m in- These cal- 
culations give lower limits, as they do not include objects 
blowing away at z > z max , but having a peak luminosity at 

Z *C Z ma x • 



4 RESULTS 

The above model allows us to determine the observational 
properties of the MIR emission from the first objects. In 
particular, we present results concerning the expected flux 



Figure 1. Expected flux as a function of the observed wavelength 
for the five H2 molecular lines considered: A e = 2.12 fim (solid 
line), 6.9 /im (dot-dashed), 9.7 /an (long dashed), 17 (short- 
dashed) and 28 fim (dotted). Fluxes are shown for halos with 
mass equal to the mean one at the relevant emission redshift (see 
text for details). 



in the various lines considered (see Table 1) and the corre- 
sponding number counts. These are then compared with the 
sensitivities of various planned instruments. 



4.1 Expected MIR Fluxes 

Fig. 1 shows the expected flux as a function of the observed 
wavelength, A , in the five different lines 2.12 fim, 6.9 fim, 
9.7 /im, 17 /im and 28 fim. As an example, we plot the fluxes 
calculated for halos of mass w 10 s — 1O 9 M0, corresponding 
to the mean one among those suffering blowaway, as given 
by the Press-Schechter formalism, at the explosion redshift. 
The curves refer to emission occurring when the shocked 
gas temperature has decreased down to values close to the 
corresponding line excitation temperature. Thus, the actual 
emission redshift is somewhat lower than the explosion red- 
shift. 

From Fig. 1 we see that the typical observed fluxes are 
in the range 10 -22 — 10~ 17 erg s _1 cm -2 , depending on the 
wavelength band (and therefore emission redshift). For the 
2.12 fim, 6.9 fim and 9.7 /jm emission lines the flux is in- 
creasing with wavelength, whereas the flux of the 17 fim 
and 28 fim excitation lines is decreasing towards higher A . 
From a general point of view this can be understood as fol- 
lows. As the IGM density increases with redshift, this causes 
a corresponding increase in the emissivity. In addition, the 
postshock gas tends to be warmer at high z, hence inducing 
preferentially the excitation of the higher excitation tem- 
perature lines. These effects are dominant for the 2.12 fim, 
6.9 fim and 9.7 fim lines, while they are overwhelmed by a 
larger value of dc and a smaller shell emission volume for 
the redder ones, producing a decreasing intensity. 

As detection chances are highest for the 2.12 fim line, 
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Figure 2. (a) Luminosity evolution for a blown-away object at 
z = 4 (long-dashed line), z = 6 (solid), z = 8 (dotted) and z. = 10 
(short-dashed). The emission wavelength is A e = 2.12 fim and the 
luminosity is evaluated for the average halo mass (see text), (b) 
Fraction of the SN mechanical energy emitted in the A e = 2.12 (im 
line as a function of time. Notation is the same as (a) . 



we have derived its luminosity evolution and the correspond- 
ing fraction of the SN energy radiated in this line by objects 
at different redshifts. As for Fig. 1, the emission is derived 
for an average halo mass. In Fig. 2a we show the luminosity 
evolution for selected redshifts in the range 4 ^ z ^ 10. As 
the explosion reshift is increased, the peak of the emission is 
shifted at earlier times and it becomes narrower. Both effects 
are due to the faster evolution of the temperature driven by 
the enhanced (i.e. oc (1 + z) 4 ) Compton cooling rates. The 
more pronounced peak occurring at higher z is due to the 
(1 + z) 3 density increase; this trend is partially counterbal- 
anced by the decreasing average halo masses with redshift. 



Figure 3. Number counts from the H2 A e = 2.12 /im emission 
wavelength as a function of the minimum observed flux, F m i n . 
The number counts are calculated for different observed wave- 
length bands with R = 3 and centered at A = 15 /im (solid 
line) and 20 fim (dashed). The limiting sensitivity of NGST in 
the same bands is also shown. 

Fig. 2b shows the time dependence of the cumulative frac- 
tion of the SN mechanical energy f e (t) = E e (t)/E emitted 
in the A e = 2.12 /im line. This quantity increases over a time 
scale inversely related to the explosion redshift and reaches a 
plateau, determined by the position of the L(i) peak seen in 
Fig. 2a. At high redshift up to 10% of the SN mechanical en- 
ergy is carried away by the considered roto- vibrational line; 
at lower redshifts this value rapidly decreases below 1%. 

4.2 Number Counts and Detectability 

By using the expressions given in eq. 12 we calculate the 
expected number of objects per unit sky area as a function 
of the limiting flux, F m i n , of a given experiment. The cor- 
responding curves are shown in Fig. 3 for the 2.12 fim H2 
line. The number counts have been estimated in two differ- 
ent MIR bands in the sensitivity range of planned imaging 
detectors on board of the Next Generation Space Telescope^ 
(NGST). For our study, we consider the R = 3 imaging mode 
of NGST, in the two bands centered at A D = 15 and 20 fim. 
In the 20 /im band, the largest fluxes are at the 6 x 10~ 18 
erg s _1 cm -2 level, where about 10 4 objects/deg 2 should be 
seen; the surface density of sources increases rapidly as F m i n 
decreases. In the 15 /tm band objects are fainter and a rapid 
drop is seen in the number counts above 10~ 18 erg s^ 1 cm" 2 . 

In the above bands, NGST will reach a sensitivity of 
~ 4 x 10 -18 erg s _1 cm -2 (shown in the Figure) for a 3a de- 
tection in one hour. At the above sensitivity level, and with 
its field of view of 4'x4', NGST should be able to detect 
about 200 sources in the 20/im band in one hour integration 

t http://augusta.stsci.edu 
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at the 3a level. These objects are all located in the redshift 
range Az ~ 7 — 10, a very intriguing time during cosmic evo- 
lution, perhaps bracketing the reionization epoch. Assuming 
a diffraction limited telescope (angular resolution ~ 0.25") 
at this wavelength, field crowding should not be a problem, 
the number of sources being well below the confusion limit. 

Longer exposure times might also allow detections in 
the 15 /jm band. Sensitivities comparable to the NGST ones 
are also expected for the SPICA HII/L2 future mission^. 
Other MIR planned facilities, such as SIRTF or FIRST will 
only reach sensitivities ~ 10~ 14 — 10~ 15 erg s _1 cm -2 , de- 
pending on observational bands. 

In conclusion, NGST appears to constitute the best tool 
to reveal the emission from the first objects; while other 
above mentioned instruments will not reach the required 
sensitivities. 



5 CONCLUSIONS 

We have investigated the possibility to detect the first ob- 
jects in the universe through the MIR line emission asso- 
ciated with the stellar feedback (i.e. blowaway of the gas) 
induced by their supernovae. A noticeable fraction of the 
explosion energy (up to 10% in the 2.2 /im line alone) is 
eventually converted and carried away by molecular hydro- 
gen roto-vibrational lines. By using a detailed treatment of 
the non-equilibrium formation and emission of H2 molecules 
behind these cosmological blastwaves we have calculated the 
observed fluxes and number counts of primordial objects 
and compared them with the expected sensitivities of future 
instruments. At the limiting flux of NGST, we do predict 
that about 200 primordial objects can be detected 20 /im in 
barely one hour of observation time. 

This detection would allow to directly image the first 
star formation in the universe and our results show that 
MIR observations, when compared with estimates of previ- 
ous studies in other bands and/or exploiting different strate- 
gies, represent a superb tool for this study. In addition, the 
proposed excitation mechanism will also allow to test and 
calibrate the stellar feedback process in the dark ages and 
finally assess the degree at which the first galaxies (and the 
IGM) have been influenced by supernova energy deposition. 
This obviously holds the key for the understanding of the 
subsequent evolution and formation of larger galactic blocks 
in most cosmological models. MIR line emission allows the 
high redshift universe (z £ 5) to be much more easily ex- 
plored than in the near IR bands, where the intergalactic 
absorption might be found to strongly affect the detectabil- 
ity of objects located close or beyond the reionization epoch: 
for example, the predicted dust opacity to sources located at 
redshift ~ 5 is as high as ~ 0.13 at the observed wavelength 
An ~ 1pm, and could considerably affect observations of the 
distant universe in that band (Ferrara et al. 1999; Loeb & 
Haiman 1997). It is worth noting that our predictions con- 
stitute lower limits to the number of observable objects as 
we have not considered possible density enhancements de- 
scending from explosions occurring inside overdense regions, 
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as the filaments of the cosmic web (and also because of some 
details of the calculations, see Sec. 3). 

Molecular hydrogen lines might also be excited during 
the collapse of the so-called PopIII objects, which rely on 
this molecule to collapse and form stars. However, there are 
several reasons to suspect that the contribution of this pro- 
cess will be negligible with respect to the emission due to the 
stellar feedback. First, if a UV background is present, molec- 
ular hydrogen can be destroyed by photodissociation. Even 
if there is no UV background, dissociation of molecular hy- 
drogen by internal UV radiation emitted from massive stars, 
formed in the high density regions of the objects, is very ef- 
ficient and the evolution of the objects and star formation 
are strongly affected (Omukai & Nishi 1999). This regulation 
effect is efficient as long as the line emission of molecular hy- 
drogen is the main cooling process, i.e. if Z ,$ 10~ 2 Z Q and 
T < 8000 K (Nishi & Tashiro 2000). Second, the mass of 
these objects is very small and therefore the amount of gas 
in the required thermodynamic state for the emission is very 
limited. Interestingly, the relative little importance of PopIII 
objects in terms of their radiation power in ionizing photons 
with respect to objects of mass above Mh has already been 
established by Ciardi et al. (2000). The latter objects emit 
a very large fraction of their binding energy in the hydro- 
gen Lya line which is therefore a much better tracer of their 
formation (Haiman, Spaans & Quataert 2000). 
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